Estimating outflow facility through pressure dependent pathways of the human eye

We develop and test a new theory for pressure dependent outflow from the eye. The theory comprises three main parameters: (i) a constant hydraulic conductivity, (ii) an exponential decay constant and (iii) a no-flow intraocular pressure, from which the total pressure dependent outflow, average outflow facilities and local outflow facilities for the whole eye may be evaluated. We use a new notation to specify precisely the meaning of model parameters and so model outputs. Drawing on a range of published data, we apply the theory to animal eyes, enucleated eyes and in vivo human eyes, and demonstrate how to evaluate model parameters. It is shown that the theory can fit high quality experimental data remarkably well. The new theory predicts that outflow facilities and total pressure dependent outflow for the whole eye are more than twice as large as estimates based on the Goldman equation and fluorometric analysis of anterior aqueous outflow. It appears likely that this discrepancy can be largely explained by pseudofacility and aqueous flow through the retinal pigmented epithelium, while any residual discrepancy may be due to pathological processes in aged eyes. The model predicts that if the hydraulic conductivity is too small, or the exponential decay constant is too large, then intraocular eye pressure may become unstable when subjected to normal circadian changes in aqueous production. The model also predicts relationships between variables that may be helpful when planning future experiments, and the model generates many novel testable hypotheses. With additional research, the analysis described here may find application in the differential diagnosis, prognosis and monitoring of glaucoma

Bone balance within a cortical BMU: Local controls of bone resorption and formation

Maintaining bone volume during bone turnover by a BMU is known as bone balance. Balance is required to maintain structural integrity of the bone and is often dysregulated in disease. Consequently, understanding how a BMU controls bone balance is of considerable interest. This paper develops a methodology for identifying potential balance controls within a single cortical BMU. The theoretical framework developed offers the possibility of a directed search for biological processes compatible with the constraints of balance control. We first derive general control constraint equations and then introduce constitutive equations to identify potential control processes that link key variables that describe the state of the BMU. The paper describes specific local bone volume balance controls that may be associated with bone resorption and bone formation. Because bone resorption and formation both involve averaging over time, short-term fluctuations in the environment are removed, leaving the control systems to manage deviations in longer-term trends back towards their desired values. The length of time for averaging is much greater for bone formation than for bone resorption, which enables more filtering of variability in the bone formation environment. Remarkably, the duration for averaging of bone formation may also grow to control deviations in long-term trends of bone formation. Providing there is sufficient bone formation capacity by osteoblasts, this leads to an extraordinarily robust control mechanism that is independent of either osteoblast number or the cellular osteoid formation rate. A complex picture begins to emerge for the control of bone volume. Different control relationships may achieve the same objective, and the ‘integration of information’ occurring within a BMU may be interpreted as different sets of BMU control systems coming to the fore as different information is supplied to the BMU, which in turn leads to different observable BMU behaviors

Estimating outflow facility parameters for the human eye using hypotensive pressure-time data

We have previously developed a new theory for pressure dependent outflow from the human eye, and tested the model using experimental data at intraocular pressures above normal eye pressures. In this paper, we use our model to analyze a hypotensive pressure-time dataset obtained following application of a Honan balloon. Here we show that the hypotensive pressure-time data can be successfully analyzed using our proposed pressure dependent outflow model. When the most uncertain initial data point is removed from the dataset, then parameter estimates are close to our previous parameter estimates, but clearly parameter estimates are very sensitive to assumptions. We further show that (i) for a measured intraocular pressure-time curve, the estimated model parameter for whole eye surface hydraulic conductivity is primarily a function of the ocular rigidity, and (ii) the estimated model parameter that controls the rate of decrease of outflow with increasing pressure is primarily a function of the convexity of the monotonic pressure-time curve. Reducing parameter uncertainty could be accomplished using new technologies to obtain higher quality datasets, and by gathering additional data to better define model parameter ranges for the normal eye. With additional research, we expect the pressure dependent outflow analysis described herein may find applications in the differential diagnosis, prognosis and monitoring of the glaucomatous eye

Cell organisation in the colonic crypt: A theoretical comparison of the pedigree and niche concepts

The intestinal mucosa is a monolayer of rapidly self-renewing epithelial cells which is not only responsible for absorption of water and nutrients into the bloodstream but also acts as a protective barrier against harmful microbes entering the body. New functional epithelial cells are produced from stem cells, and their proliferating progeny. These stem cells are found within millions of crypts (tubular pits) spaced along the intestinal tract. The entire intestinal epithelium is replaced every 2–3 days in mice (3–5 days in humans) and hence cell production, differentiation, migration and turnover need to be tightly regulated. Malfunctions in this regulation are strongly linked to inflammatory bowel diseases and to the formation of adenomas and ultimately cancerous tumours. Despite a great deal of biological experimentation and observation, precisely how colonic crypts are regulated to produce mature colonocytes remains unclear. To assist in understanding how cell organisation in crypts is achieved, two very different conceptual models of cell behaviour are developed here, referred to as the ‘pedigree’ and the ‘niche’ models. The pedigree model proposes that crypt cells are largely preprogrammed and receive minimal prompting from the environment as they move through a routine of cell differentiation and proliferation to become mature colonocytes. The niche model proposes that crypt cells are primarily influenced by the local microenvironments along the crypt, and that predetermined cell behaviour plays a negligible role in their development. In this paper we present a computational model of colonic crypts in the mouse, which enables a comparison of the quality and controllability of mature coloncyte production by crypts operating under these two contrasting conceptual models of crypt regulation

Analysis of Wnt signaling β-catenin spatial dynamics in HEK293T cells

Background Wnt/β-catenin signaling is involved in different stages of mammalian development and implicated in various cancers (e.g. colorectal cancer). Recent experimental and computational studies have revealed characteristics of the pathway, however a cell-specific spatial perspective is lacking. In this study, a novel 3D confocal quantitation protocol is developed to acquire spatial (two cellular compartments: nucleus and cytosol-membrane) and temporal quantitative data on target protein (e.g. β-catenin) concentrations in Human Epithelial Kidney cells (HEK293T) during perturbation (with either cycloheximide or Wnt3A). Computational models of the Wnt pathway are constructed and interrogated based on this data. Results A single compartment Wnt pathway model is compared with a simple β-catenin two compartment model to investigate Wnt3A signaling in HEK293T cells. When protein synthesis is inhibited, β-catenin decreases at the same rate in both cellular compartments, suggesting diffusional transport is fast compared to β-catenin degradation in the cytosol. With Wnt3A stimulation, the total amount of β-catenin rises throughout the cell, however the increase is initially (~first hour) faster in the nuclear compartment. While both models were able to reproduce the whole cell changes in β-catenin, only the compartment model reproduced the Wnt3A induced changes in β-catenin distribution and it was also the best fit for the data obtained when active transport was included alongside passive diffusion transport. Conclusions This integrated 3D quantitation imaging protocol and computational modeling approach allowed cell-specific compartment models of the signaling pathways to be constructed and analyzed. The Wnt models constructed in this study are the first for HEK293T and have suggested potential roles of inter-compartment transport to the dynamics of signaling

Time evolution of deformation in a human cartilage under cyclic loading

Recent imaging has revealed that in vivo contact deformations of human knee cartilage under physiological loadings are surprisingly large—typically on the order of 10%, but up to 20 or 30% of tibiofemora cartilage thickness depending on loading conditions. In this paper we develop a biphasic, large deformation, non-linear poroelastic model of cartilage that can accurately represent the time dependence and magnitude of cyclic cartilage deformations in vivo. The model takes into account cartilage tension–compression nonlinearity and a new constitutive relation in which the compressive stiffness and hydraulic permeability of the cartilage adjusts in response to the strain-dependent aggrecan concentration. The model predictions are validated using experimental test results on osteochondral plugs obtained from human cadavers. We find that model parameters can be optimised to give an excellent fit to the experimental data. Using typical hydraulic conductivity and stiffness parameters for healthy cartilage, we find that the experimentally observed transient and steady state tissue deformations under cyclic loading and unloading can be reproduced by the model. Steady state tissue deformations are shown to cycle between 10% (exudation strain) and 20% (total strain) in response to the cyclic test loads. At steady-state cyclic loading, the pore fluid exuded from the tissue is exactly equal to the pore fluid imbibed by the tissue during each load cycle

Colon cryptogenesis: Asymmetric budding

The process of crypt formation and the roles of Wnt and cell-cell adhesion signaling in cryptogenesis are not well described; but are important to the understanding of both normal and cancer colon crypt biology. A quantitative 3D-microscopy and image analysis technique is used to study the frequency, morphology and molecular topography associated with crypt formation. Measurements along the colon reveal the details of crypt formation and some key underlying biochemical signals regulating normal colon biology. Our measurements revealed an asymmetrical crypt budding process, contrary to the previously reported symmetrical fission of crypts. 3D immunofluorescence analyses reveals heterogeneity in the subcellular distribution of E-cadherin and β-catenin in distinct crypt populations. This heterogeneity was also found in asymmetrical budding crypts. Singular crypt formation (i.e. no multiple new crypts forming from one parent crypt) were observed in crypts isolated from the normal colon mucosa, suggestive of a singular constraint mechanism to prevent aberrant crypt production. The technique presented improves our understanding of cryptogenesis and suggests that excess colon crypt formation occurs when Wnt signaling is perturbed (e.g. by truncation of adenomatous polyposis coli, APC protein) in most colon cancers

Theoretical analysis of the spatio-temporal structure of bone multicellular units

Bone multicellular units (BMUs) maintain the viability of the skeletal tissue by coordinating locally the sequence of bone resorption and bone formation performed by cells of the osteoclastic and osteoblastic lineage. Understanding the emergence and the net bone balance of such structured microsystems out of the complex network of biochemical interactions between bone cells is fundamental for many bone-related diseases and the evaluation of fracture risk. Based on current experimental knowledge, we propose a spatio-temporal continuum model describing the interactions of osteoblastic and osteoclastic cells. We show that this model admits travelling-wave-like solutions with well-confined cell profiles upon specifying external conditions mimicking the environment encountered in cortical bone remodelling. The shapes of the various cell concentration profiles within this travelling structure are intrinsically linked to the parameters of the model such as differentiation, proliferation, and apoptosis rates of bone cells. The internal structure of BMUs is reproduced, allowing for experimental calibration. The spatial distribution of the key regulatory factors can also be exhibited, which in diseased states could give hints as to the biochemical agent most accountable for the disorder

Yield stress measurements of aqueous foams in the dry limit

This paper reports measurements of yield stress of aqueous foams approaching the dry foam limit using a pendulum device. Traditionally, the vane rheometer has been used to measure the yield stress in liquids that exhibit wall slip. However, using the simple and inexpensive pendulum technique, shear rates many orders of magnitudes lower can be achieved. The pendulum was used to observe the change in yield stress for the foam as the gas fraction and bubble size increased. The local gas fraction in the foam was found by measuring the sonic velocity, and the bubble size was determined photographically. Strong support is found for the existence of a true yield stress in aqueous foams at the dry foam limit. Yield stress results, once scaled by σ/⟨R⟩, agree well with data from previous studies

Micro-computed tomographic analysis of the radial geometry of intrarenal artery-vein pairs in rats and rabbits: Comparison with light microscopy

We assessed the utility of synchrotron-radiation micro-computed tomography (micro-CT) for quantification of the radial geometry of the renal cortical vasculature. The kidneys of nine rats and six rabbits were perfusion fixed and the renal circulation filled with Microfil. In order to assess shrinkage of Microfil, rat kidneys were imaged at the Australian Synchrotron immediately upon tissue preparation and then post fixed in paraformaldehyde and reimaged 24 hours later. The Microfil shrank only 2-5% over the 24 hour period. All subsequent micro-CT imaging was completed within 24 hours of sample preparation. After micro-CT imaging, the kidneys were processed for histological analysis. In both rat and rabbit kidneys, vascular structures identified in histological sections could be identified in two-dimensional (2D) micro-CT images from the original kidney. Vascular morphology was similar in the two sets of images. Radial geometry quantified by manual analysis of 2D images from micro-CT was consistent with corresponding data generated by light microscopy. However, due to limited spatial resolution when imaging a whole organ using contrast-enhanced micro-CT, only arteries ≥100 and ≥60 μm in diameter, for the rat and rabbit respectively, could be assessed. We conclude that it is feasible and valid to use micro-CT to quantify vascular geometry of the renal cortical circulation in both the rat and rabbit. However, a combination of light microscopic and micro-CT approaches are required to evaluate the spatial relationships between intrarenal arteries and veins over an extensive range of vessel size